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Abstract 


This report describes a system of computer programs for testing 
Iterative methods for solving nonlinear systems of equations. The 
programs are written In FORTRAN V for the Unlvac 1108. The system can 
be updated by adding or deleting methods and sets of equations, and It 
can be used with various sets of Initial approximations for the Iterations. 
The results of the Iterations are collected In a data file which can be 
printed in table form for comparing different methods. 
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1 . Introd uct ion 

During the pest decade growing attention ha* been directed toward *he 
development of Iteiatlve methods for the numerical solution of systems of 
nonlinear equations 

i 

f ] (x 1 x n ) = 0 

f 2 ( x j • • • • , x n ) = 0 

(D 

f n (Xi , . . . ,x n ) = 0 

For example, besides the many variations of Newton's method, there are 
several variations of the secant method, the so-called generalized linear 
methods, and the large number of minimization methods in case the f.'s are 
the partial derivatives of a nonlinear functional. For a survey of these 
methods we refer to the forthcoming book by J. M. Orteqa and W. C. Rhelnboldt*. 

None of these methods Is satisfactory for all types of systems, and, in 
fact, for different systems the behavior of a certain method may change dras- 
tically. In addition, the convergence behavior of any method depends, In most 
cases very critically, on the choice of the Initial approxlmation(s) . Although 
there Is a growing literature about convergence theorems for the different 
methods, the conditions of many of these theorems for a specific system are 
frequently not numerically checkable or they predict only very small conver- 
gence domains. Moreover, only a few results are known about the stability 
behavior of the various methods under round-off. 

^ Iterative Solution of n-Dimensional Nonlinear Equations , Blalsdell Publishing 
Ci 
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More than In many other areas of numerical analysis, it Is therefore 
of considerable Importance to conduct computational experiments with 
various methods, in each case, for many different choices of systems and 
for a large number of initial approximations. A series of such experiments 
with one given method may give sorne Insight Into the behavior of that method 
under different conditions. Usually, It is more Important to have an exper- 
imental comparison of different methods under varying conditions. 

Such a test program will quickly run into considerable data-handl 1 ng 
difficulties unless from the very outset a flexible Integrated system of 
programs has been devised which allows the Incorporation of new methods, 
systems and Initial approximations at any time, and which permits the 
accumulation of test results In a form suitable for future processing. 

This report describes suen a system of programs for testing any number 
of iterative methods fo<" solving any number of systems of equations (1) 
under arbitrary Initial conditions. In the next chapter, we describe the 
adopted numerical tests for the convergence behavior of an iterative 
method. The third chapter describes the various files of the system, and 
the fourth chapter, the test program. The last two chapters give details 
of the presently available postprocessing programs. 
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2 . Numer ica l tests for c onver gence 
Let 


( 2 ) 


X f X , . . . j X , X 


>e the sequence of Iterates obtained from some process after i steps, and 


(3) 


Fq 7^1 .... *Fi_i ,Fi 


tl.s corresponding function values. Here 

= (x:j,x^,. .. ,xj*) T , Fj * (t' 1 (x^),t 2 (x'^),...,f n (x^)) T 

are n-d-imenslonal vectors with the norms II x^H and || F ^ li . 

The convergence tests are to be Incorporated in each different method 
program and consist of tests performed after every step of the Iteration 
to determine whether the process should terminate, and of tests performed 
on the last-obtained iterate after termination. 

The tests for the termination of an iterative process after the i-th 
step are as follows: 

a) If the method breaks down during the 1-th step, e.g., If a singular 
system of equations arises, then the Iteration is terminated and the 
termination symbol B Is generated. 

b) If iix 1 HrMO 20 or || F^ || > 10 20 , the iteration Is terminated. The 
iteration is considered to be diverging, and the termination symbol 
D Is used. 

c) If lx 1 -* 1 ’ 1 11 -^ 2 ’ the 1terati on is terminated; It is considered to 
be converging, and the termination symbol C Is generated. 

d) If 
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||x , -x i ' 1 IMIx i ‘ 1 -x 1 ' 2 l| < ...<||x 1 * 1 ® 4l -x 1 ‘ 1 °ll 

and (1 > 1 Q ) 

llx^x 1 - 1 !! < Ux 1 !! , 

then the iteration is terminated; it is considered to be converging, 
and again the symbol C is generated. 

e) If 

#x 1 -x 1 - , ll>llx i - , -x , - 2 ll> ... Ix 1 ' 1 0 +1 -x 1 ' 1> | , (1 > 1 0 ) 

then the iteration is terminated; it is considered to be diverging; 
this is again denoted by tf. 

f) If 

IlF. IJ> II F . _ 1 H > ... > Hf._. +1 H 

and (1 >. 1 Q ) 

llx 1 - X 1 ' 1 ll ^ !lx 1 " 1 -x 1 “ 2 1) , 


then the iteration is terminated; it is considered to be diverging 
and also indicated by D. 

g) If i > MAX, the iteration is terminated; It is considered to be indecisive, 
and the termination symbol I is generated. 

The second group of tests a*"e performed when the termination symbol was 
either B (case a) or C (case c and d). If 

II r i || < € 1 
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then B is changed to BC. If 


II F t || >£ 1 

then C is changed to CB. BC indicates an instability of the method 
near the solution. CB indicates convergence to a point which is not a 
solution of (1) or very slow convergence. 

The parameters, MAX, 6^, 6 2 , and Iq* a ^ e input parameters; 
they can be changed for different iterations. 
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3 . Fil e organization 

The Meth od-file is the collection of various iterative methods for 
solving any nonlinear sys .em of equations. Each method has a unique n. 
consisting of not more than six characters. 

The Functi on-file (or F-f lie ) is the collection of nonlinear systems 
of equations ef the form (1). These systems of equations are sequentially 
numbered in the file, and one specific system of equations is specified 
by the letter F and its sequence number. 

The Initial-value-file (or S-f lie ) is the collection of sets of initial 
approximations for the iterative processes. These sets are sequentially 
numbered in the file, and k-S refers to the k-th set of initial approxi- 
mations. Each set of initial approximations consists of one or several points 
in R n . These are again numbered sequentially, and k-S-i refers to the 
i-th point In the set k-S. 

The PR-file is the program file. All programs associated w'th the 
Test System were written in FORTRAN V* for the Uni vac 1108. The symbolic 
and relocatable forms of these programs are collected on this t?pe-file 
including the programs associated with the previously-defined files. The 
file may be updated by adding or deleting method subroutines, function rou- 
tines, initial-value generating subroutines, cr with any auxiliary programs. 
Using the collector of EXEC 8 for the Univac 1108, any subprogram can be 
selected from the PR-file for execution. The following rules for naming 
programs must be applied when the file is updated: 



♦Except a small assembly language routine which gives the elapsed time 
msec. 
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Symbolic 

Relocatable 

Method subroutine 

A. METHOD/' NAME' 

A.METHDX/'NAME 

Initial value generating 
subroutine 

A.S/n 

A.SX/n 

Function routine 

A.F/m 

A.FX/m 


A. D/m 

A.DX/m 


where 'NAME' is the unique name of the method, and n and m are the appro- 
priate sequence numbers. Further description of these routines can be found in 
4.5 - 4.7. 

At the present, the PR-file contains three main programs: 

i. Method Testing Program 

ii. List of Tape Output File 

iii. Table Generating Program 

The Method Testing Program, consisting of several subprograms, tests one 
method Dy solving one system of nonlinear equations with one set of initial 
approximations. The List of Tape Output File lists the contents of File 2 
(see below). The Table Generating Program collects data from File 2 and prints 
it out in table form. 

The Outpi't-File is the collection of the outputs of the Method Testing 
Program. It consists of printed, punched and tape outputs. The generated 
and updated tape output is called File-2. 
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4 . Meth od T esting Program 
4.1 Deck s etup 

For the execution of the Method Testing Pronrnm the following card 
deck is necessary: 


Initial 

Control 

Cards 


'RUN 

'ASG.T PR,8C,tl 
'ASG.T 2,8C,t2 

‘ASG.T A.F 

'COPIN, SR PR. 

'FREE PR. 


Cards to 
define 
one run 


i 




'MAP, I A.PROC.A.PROCX 

SEG PRO 

IN A.MAINX,A.PCHX,A.PRNTX,A GAUSX.A.SCALX.A.TIMEX 

IN A.LTWOX 

IN A.METHDX/'name* 

IN A.FX/n,A.DX/n 

IN A.SX/m 

' XQT A.PROCX 

(one data card) 


Terminating, 

cards 


'FREE 

'FREE 

'FIN 


2 . 

A. 


Initial Control Cards: These cards define the files; tl and t2 should be 

replaced by the tape reel numbers of the files PR and 2. 

The "Cards to define one run" may be repeated before the "Terminating Cards" 
to execute more than one run with different combi lations of methods, 
functions, etc. The 'IN' cards define the selection of the necessary 

programs. The first 'IN' card defines the main program and subroutines 

needed for any run. The second 'IN' card defines the subroutine which 
calculates tna norm of a vector. The third 'IN' card selects the method 
used for this run; 'name' should be replaced by the alpha-numeric name 
of the method. The fourth and fifth 'IN’ c^rds select the system of 
equations and the initial value sets; ' n ' and 'm' sho'^d be replaced by 
the corresponding file indices. 
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The "Data Card" has the following format: 


col . 1-5 

10 


14 - 15 


(Format: 


16 - 25 
26 - 35 
36 - 45 
46 - 55 
315, 4E10.5) 


Example: 

50 0 


MAX = maximum number of iterations allowed 
NPR = print index 

step-by-step print during iteration 
only final results to be printed 
NPU = Punch and Tape index 
NPU = 0 or 2 no punch 

or 3 punch final results 
0 cr 1 no tape output 
2 or 3 final results written on File 2 
negative sign indicates a new tape; 
positive sign indicates that the tape 
must be positioned. 




NPR = 1 
= 0 


= 1 


input parameters described 
in Chapter 2 


-2 


.IE-6 


1 E-6 


. 1 E-5 


. 5E1 


4.2 Available subroutines 

The following subroutines are incorporated in the main program 
package; they may be used in any new method subroutine with the following 
calling sequences: 

VMPV (N,A,IA,B,IB) is a function routine for computing the scalar product of 
the N-dimensional vectors A and B. IA and IB are the index incre- 
ments for calculating the locations of the components of A and B. 

w 

VMPV ( N , A , I A , B , I B ) = A(l+(i-l )*IA)*B(l+(i-i )*IB) 

The routine calculates the products and the summation in double precision 
arithmetic and truncates the result to single precision. 


- 10 - 


GINV(*,A,B,N) is a matrix inversion subroutine. The N*N matrix 
A is inverted by Gauss elimination with pivoting and the inverse is placed 
in B. Both matrices A and B should be dimensioned 25* 25. The input 
matrix A remains unchanged unless B refers to the same storage, which 
is permitted. The internal calculation is done in double precision arith- 
metic; * is the error return when the matrix A has been found numerically 
singular. 

GSLV (*,E,F,N) is a subroutine for solving the linear system of equations 


t 


E ( I , J )*X (J ) = F ( I ) 


The solution X replaces F and the matrix E remains unchanged; E must 
be dimensioned 2bx25. The routine solves the equations by Gauss elimina- 
tion with pivoting in double precision arithmetic and gives truncated 
final results. The return to * occurs if the matrix E is numerically 
singular. 


4.3 Outputs 

4.3.1 - Printed output 

The first page of the output is a heading page. It consists of the data 
supplied by the data card, and, of the initial value and function file identi- 
fications. An example of this heading is given on page A-l . 

Following the heading page, each iteration processed with an "Initial- 


Value" is printed as follows: 
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(a) The "Initial Value" with its index numbers in the given set S/m 

(b) Step-by-step results of the iteration consisting of the step- 
number, the maximum norm of the remainder, and the independent 
variable. If the print index NPR is zero on the data card, this 
printing is suppressed. 

(c) Results of the iteration and the file identifications. The 

results of the iteration consist of the termination symbol, the 

solution index, the number of function evaluations, and the time 

in msec spent in the method program. As described in Chapter 2, 

the termination return can be one of the following: 

C = Convergence return from the method subroutine with X as 
a last iterate and either 



or 

IIf(x) II £6, 

where Z. is one of the solutions supplied by the function 
routine 1 and the norm is defined by the selected norm 
routine. In the first case the solution index is set to 
i; if the first condition is not satisfied while the second 
one is, then the solution index is set to zero. 

CB = Convergence return from the method subroutine with X 
as a last iterate but wit.h neither one of the above conditions 
satisfied. 

D = Divergence return from the method subroutine. 

I = After the allowed maximum number of iterations neither 
convergence nor divergence was established by the method 
subroutine 

B = The method broke down and the last obtained point X did 
not satisfy either one of the conditions under 'C'. 
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BC = The method broke down but the last obtained point X 
satisfied one or both of the conditions under 'C'. 

The solution index is set to i / 0 whenever the last obtained 
point is close to the i-th solution supplied by the function 
routine, that Is, when the first condition of 'C' is satisfied. 
Otherwise, it is always s»t to zero. 


An example for the front page is shown on page A-l ; examples for the itera- 
tions are shown on pages A -2 and A-3. 


4.3.2 Punched output 

When punch is requested by the data card, one summary card is punched 
for each iteration with the following format: 


1. 

col 

. 2 




Identifier number = 4 

2. 

II 

4 

- 

.9 

= 

Name of the method 

3. 

II 

10 

- 

14 

= 

File index of the Function routine 


II 


15 


= 

Letter F 

4. 

II 

16 

- 

20 

= 

File index of the Initial Value Generating routine; m 


II 


LI 


= 

Letter S 

5. 

II 

22 

- 

24 

= 

Initial value index in the set S/m 

6. 

II 

27 

- 

28 

= 

Type of return 

7. 

II 

29 

- 

31 

= 

Solution index 

8. 

II 

32 

- 

37 

= 

Number of executed iteration steps 

9. 

II 

38 

- 

43 

= 

Number of function evaluations 

10. 

II 

44 

- 

49 

s 

Time spent in msec 

11. 

II 

51 

- 

56 

= 

Norm used 

12. 

I* 

57 

_ 

60 

= 

Maximum number of iteration steps allowed 
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13. " 61 - 66 = £ ] 

14. " 67 - 72 = £ ^ 

16. " 73 - 78 = E 3 

16. " 79 - 80 = i Q 

Page A-4 shows an example 


^ in E6.1 formats 
ing a list of output cards. 


4.3.3 Tape outpu t (File 2) 

When tape output is requested by the data card, one logical binary 
record of 26 words is written on tape (File 2) for each iteration. The 
fir:t 16 words are the 16 data listed in the Punched output; the last 10 
words are reserved for possible other classifications or remarks. A 
dummy-record, with 26 words of zeros, is use', to mark the end of the 
file, 


4.4 Norm-Function Routines 

The Method Testing System provides two norm routines: 

* ? 1/2 

L-TWO: Dxli = ( £ x. ) 

i* I 

L-MAX: II x II = Max llx.l , i=l N) 

but it is possible to insert other norm-function routines. A 
properly programmed Norm routine is illustrated on page A-5. The first 
entry FNORM (X,N), (line 2), is used to return the value of the norm of the 
N-dimensional vector X. The second entry FNORMM (FF), (line 16), is used 
to return the alphabetic name of the norm; FF is a dummy argument. 
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4.5 Method Routines 


The Testing Program requires the Method routines to be programmed 
in a standard form. This form is illustrated by the Newton Method on 
page A -6 : 

Lines 1 - 14, initialization part, and lines 28 - 49, testing part, 
are standard. 

The first line is a Fortran control statement. Note that the name 

of the method appears as a version. 

The arguments of the routine are as follows (line 2): 

*1 Return position if divergence is detected 

*2 Return position if no divergence or convergence after 

MAX number of iteration steps 

*3 Return position if the method breaks down. In the 

example, this return occurs when the Newton method 
tries to solve a numerically singular linear system 
(lines 22 and 52) 

N ( Dimension 

X'.i) i = 1,.. ,N. Initial value vector, which is replaced by 
the i cerated values 

MAX Maximum number of iterative steps allowed 

ACC ( 1 ) £,. "zero" for the norm of the dependent variable, F (X ) 

ACC (2) ei. "zero" for the norm of the independent variable, X 

ACC '3) 6 3 . relative accuracy required for X 

ACC(4) ip, in floating point format, number of consecutive terms 
required to detect convergence or divergence 
IT Number of iterations performed 

AN Name of the method 

In linps 8-11, the counters are set to zero; IT - iterative step, 

ICC = convergence counter, IDC1 = divergence counter for X, IDC2 = divergence 
counter for F(X). 

In lines 12 - 14, i Q as an integer is established, the output routine 
is called with the initial value and the name of the method is placed into 
the argument. 


Lines 15 - 27 are dependent upon the particular method, except fo»* 
line 20, where the iteration loop starts. 

The testing part (line? 28 - 49) assumes the newly-obtulned Independent 
variable Xj T In X and the last correction (X IT - X j j_ -j ) or (aj T _-| - X JT ) 
in DX. First the Output routine is called, then the norms are calculated. 

If the norm of the last correction is less than a normal (C) return 

occurs. To avoid machine overflow, a divergence return occurs if the norm 
of X or F (X ) is greater than or equal to 10^° (line 33). Lines 35 - 40 
update the counters except at the first iteration step. Lines 42 - 46 test 
for convergence or divergence; they are skipped in the first ij» steps. 

Lines 47 and 48 save the norms for the next iteration step, and line 49 tests 
the total number of iterations. 

The special comment lines (16, 27, 50) indicate where statements 
particular to the method may be inserted. 

4.6 Function Routines 

A properly programmed function routine is illustrated on page A-7. 

The Fortran control lines, 1 and 19, contain the File index number of 
the function as version (1 in this example). 

The first entry FXK(K,X), line 2, is used to calculate the remainder 
of the k-th equation at the given point X, and the counter NF for the 
number of function evaluations Is incremented. 

The second entry, NFEK(K) , line 13, gives the number of function 
evaluations and resets the r ounter; K is a dummy argument. NF and FN 
are made equivalent to avoid floating point-integer conversion. 
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The third entry, FX(F,X), line 20, gives the remainders of all equations 
at the given point X. 

The fourth entry, DERI VE (X ,DF) , line 26, gives the Jacobian. DF ; 
at iiie given point X and increments NDV by l? for the number of 
function evaluations. N is the dimension. 

The fifth entry, NFCV(NE), line 36, sets NE to the final number of 
function evaluations and resets the counters, NF and NDV. 

The sixth entry, FHEAD(AL,INDF N,NX,XX), line 42, gives the necessary 
information about the system of equations. INDF is set to the file index 
number, N is set to the dimension, the 26 words of AL are set to a 
heading, NX is set to the number of solutions given in XX, and finally, 
XX(I,J); I = 1,N; J = 1,NX is set to the solutions. 

4_._7 I nitial-Value Generating Routine 

A properly programmed Initial-Value Generating Routine is illustrated 
on page A-8. The Fortran control statement (line 1) uses the file index 
number (1 in this example) as a version. The first entry INVAL(*,X), line 2, 
is used to return a new n i t i a 1 value in X upon subsequent calling. A 
return to * occurs when all initial values generated by the routine were 
given before. The second entry INVIND(INDS) , line 26, is used to return 
the file index number in INDS. 

The presented example is a type of Initial-Value Generating Routine 
which we adopted for 2-dimensional cases. The generated initial values 
are on concentric circles. The origin of the circles is returned by using 


- 17 - 


the third entry (PARM(M,Z), line 29. The initial value points on one circle 
are equally distributed, but the circles can be rotated so that the correspond- 
ing points on the different circles will lie on a spiral. The parameters 
are listed in a Data Statement (line 11). By changing this statement, it 
is easy to obtain another routine for the 2-dimensional case; naturally, 
the new routine should have a different file index. 
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5. List of Tape-Output File Program 

This r 'rogram lists the contents of the Tape-Output File (File 2). 

For the outputs of each run by the Method Testing Program, the List 
program prints one line of output. An example is shown on page A-9, 

The notation for the Initial-Value File is as follows: the first integer 

is the index number of the Initial -Value Generating Program. It is 
followed by the letter S. The next two integers indicate the numbering 
of the initial values in the set generated by the routine, i.e., 1 SI -40 , 
indicates the first Initial-Value Generating program which generated forty 
different initial values numbered from 1 to 40. 


Deck Set-up: 


'RUN 

f 0 • 

'ASG.T 

PR,8C,tl 

‘ASG,T 

2,8C,t2 

'ASG J 

A,F 

•COPIN 

PR, A 

•FREE 

PR 

•XQT 

•FIN 

A.LISTX 


6. Table Generating Program 


The Table Generating Program selects data from the Tape Output file 
and arranges them in a table form. This printed table provides a quick, 
visual comparison of different methods. A sample output is shown on page 
A-10. 

There are three types of selections to be defined by two data cards: 

1. The basic selection defines, for all iterations to be included in 
the table, the indices of the function file and of the initial-value file, 
the maximum allowable number of iterations, the norm, and the number i^ 
of consecutive terms used in the convergence and divergence criteria. 

2. The method selection defines the methods by name. A maximum of 
nine methods can be used for one table, and they define the columns of the 
table. The rows of the table are defined by the different initial values 
in the set of initial values defined by the basic selection. 

3. The result selection defines the data to be printed in the table. 

These always consist of the type of return and nothing else or one of the 
following: the number of iterations, the number of function evaluations, 

the time spent in msec, and the solution index. 

The format of the two data cards is as follows: 

Function file index 
Initial Value file index 
Maximum number of iterations 
Name of the norm 


First card: col . 1 - 5 

" 6-10 
" 11 - 15 

" 20 - 25 

■ 31 - 35 
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Second card: col. 1 - 5 

" 8-13 

" 16 - 21 

" 24 - 2S 


Number of methods to be listed 
Name of the first method 
" " " second " 

" " " third " 


up to 
" 72 - 77 


ninth " 


where ND is defined as the result selection: 


if NO = 1 only the type of return is to be printed 


2 type of return and the number o* iteration 



2 ii ii H 

" the 

solution index 


^ it ii ti 

" the 

number of function evaluations 


5 

" the 

time spent in msec 

Deck Set-up: 




'RUN 




1 ASG,T 

PR*8C,tl 



1 ASG,T 

2 ,8C,t2 



'ASG,T 

A,F 



'COPIN 

PR, A 



'FREE 

PR 



' XQT 

A.TBX 



(set of 2 data cards) 




'FIN 


MAX. NuMfaER CF ITERATIONS 

Z 

50 

Print iNCtx 

— 

1 

Punch INDEX 

— 

3 

ACCUKACY REQUIREMENTS 

— 

.1—06 *1— Ob .1—05 .5+01 

norm 

z 

L-TWO 

INITIAL VALUE FILE 

z 

IS 

Dimension 

— 

2 

function file 

• . ^ i r— » a \ 

- 

IF 


( IF ) F ( 1 ) r x ( 1 ) **2 + X ( 2 ) **2 - l.o 
F ( 2 ) = X(D**3 - X<2) - 1.0 


3 KNOWN SOLUTIONS = 

1 .iooooouo+ni .oooooooo 

2 .oooooouo lcooonoo+oi 

3 • 54368°U1— 00 -.83928675-00 








A-1 


]-TH INITIA L VALUES 

.iioouooo+oi .00000000 


0 .39-00 

1 .31-01 

2 .20-03 

3 .17-06 

4 .54-07 

5 .54-07 


• 11000000 + 01 
•10045455+01 
•10001352+01 
. 10000001+01 
. 10000000 + 01 
• 10000000+01 


.uooonooo 

. 15500005-nl 

.34720078-03 

.15390106-06 

.40894968-08 

.48894962-08 


RETURN *1TH C SOLUTION = 1 

NUMBER OF FUNCTION EVALUATION = 32 TiME = 32 vS. 


METHOD nL*TON» FUNCTION 1» IN.VAL. 1. NORM L-T*0 


A-2 



15-TH 

INITIAL VALUES 



.12878553+01 

-.52654954-00 


0 

.19+01 

•12876553+01 

-.52654954-pO 

1 

.35-00 

.98180611-00 

- . 38634059-00 

2 

.19+01 

•4505U039-00 

-.15900403-01 

3 

.39-00 

.29373206-01 

-.11649747+01 

4 
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1S1- 

40 

L-TWO 

50 

.1-06 

.1-06 

.1-05 

5 

b 

SECNTP 

IF 

1S1- 

40 

L-TWO 

50 

.1-06 

.1-06 

.1-05 

5 

6 

A3l5 

IF 

1S1- 

40 

L-TWO 

50 

.1-06 

.1-06 

.1-05 

5 

7 

A3l6 

IF 

1S1- 

40 

L-TWO 

50 

.1-06 

.1-06 

.1-05 

5 
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table 

NU^BFR O c ITERATIONS - TyPF OF REtURN 


FUNCTION FILE 


- 

i 

initial value 

FILF 

= 

i 

MAX. NUMBER OF 

ITER 

• * 

50 

NORM 


- 

L-TWO 

10 



5 


T N.VAL. 

NEWTON 

moDnwT 

secant 

secnti 

SECNTP 

AU5 

1 

5-C 

6-C 

6-C 

6-C 

6-C 

5-C 

2 

4-C 

7-C 

6-C 

6-C 

6-C 

4-C 

3 

4-C 

9-C 

7-C 

7-C 

7-C 

4-C 

4 

4-C 

7-C 

6-C 

6-C 

6-C 

4-C 

5 

5-C 

7-C 

5-C 

5-C 

5-C 

5-C 

6 

4-C 

8-C 

5-C 

5-C 

5-C 

4-C 

7 

5-C 

1 5-C 

4-C 

4-C 

5-C 

6-C 

8 

4-C 

10-C 

6-C 

6-C 

6-C 

4-C 

9 

6-C 

24-C8 

1 4-C 

1 4-C 

14-C 

6-C 

10 

6-C 

24-C 

8-C 

8-C 

8-C 

6-C 

1 1 

6-C 

50-1 

9-C 

9-C 

8-C 

6— f 

12 

8-C 

6-D 

9-C 

9-C 

12-C 

5-C 

13 

8-C 

16-0 

1 5— C 

15-C 

in-c 

6-C 

14 

6-C 

9-0 

U-C 

U-C 

9-C 

6-C 

15 

9-C 

50-1 

14— C 

1 4— C 

12-C 

8-C 

16 

19-C 

4—0 

9-0 

9-0 

35-C 

8-C 

17 

6-C 

35-CB 

12-C 

12-C 

11-C 

6-C 

18 

6-C 

45-CB 

9-C 

9-C 

9-C 

6-C 

19 

13-C 

3-D 

21-D 

21-D 

50-1 

6-C 

20 

24-C 

3-D 

6-CR 

8-C8 

23-RC 

1 7-C 

21 

7-C 

5-D 

8-C 

8-C 

7-C 

b-C 

22 

7-C 

50-1 

29-C 

27-C 

12-C 

7-C 

23 

6-C 

6-0 

12-C 

12-C 

12-C 

6-C 

24 

9-C 

6-0 

1 4— C 

1 4— C 

13-C 

9-C 

25 

7-C 

50-1 

10-C 

10-C 

12-C 

7-C 

26 

38-C 

2-0 

44-0 

44-0 

50-T 

b-0 

27 

lfl-C 

5-D 

17-C 

17-C 

27-BC 

1 7— C 

28 

13-C 

50-1 

18-C 

13-C 

18-CB 

l 3— C 

29 

6-C 

50-1 

8-C 

8-C 

9-C 

6-C 

30 

6-C 

10-D 

U-C 

U-C 

11-C 

6-C 

31 

17-C 

3-D 

7-CB 

7-CR 

36— C 

9-C 

32 

9-C 

1 1-D 

1 4— C 

14-C 

14-C 

9-C 

33 

18-C 

3-0 

14-0 

1 4— D 

24-C 

8— C 

54 

U-C 

6-D 

46— CR 

47— CB 

12-C 

6-D 

35 

16-C 

50-1 

50-1 

24-BC 

29-RC 

1 6-C 

36 

20-C 

3-0 

2 7-C 

24-C 

28-C 

17-C 

37 

8-C 

50-1 

20-C 

18-C 

27-C 

8-C 

38 

8-C 

6-D 

30-8C 

27-8C 

23-C 

8-C 

39 

7-C 

6-0 

17-C 

17-C 

17-C 

7-C 

40 

7-C 

45-CR 

12-C 

12-C 

12-C 

7-C 
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